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Abstract The narrowness of the distribution of the peak energy of vF v 
spectrum of gamma-ray bursts (GRBs) and the unification of GRB pop- 
ulation are great puzzles yet to be solved. We investigate the two puzzles 
based on the global spectral behaviors of different GRB population, the long 
GRBs, the short GRBs, and the X-ray flashes (XRFs), in the HR — E P plane 
(HR the spectral hardness ratio) with BATSE and HETE-2 observations. 
It is found that the long GRBs and the XRFs observed by HETE-2 seem 
to follow the same sequence in the HR — E p plane, with the XRFs at the 
low end of this sequence. We fit the sequence by a universal Band function, 
and find that this sequence is mainly defined by the low energy index a, 
and is insensitive to the high energy index, (3. With fixed (3 = —5, a best 
fit is given by a — —1.00 with Xmin/dof = 2.2. The long and short GRBs 
observed by BATSE follow significantly different sequences in the HR — E p 
plane, with most of the short GRBs having a larger hardness ratio than the 
long GRBs at a given E p . For the long GRBs, a best-fit gives a = —0.30 and 
[3 = -2.05. For the short GRBs, a best fit gives a = -0.60 with xLn = I- 1 
(with (3 fixed at —2.0 because it is numerically unstable). The a value for 
the short GRBs is significantly greater than that for the long GRBs. These 
results indicate that the global spectral behaviors of the long GRB sample 
and the XRF sample are similar, while that of short GRBs is different. The 
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short GRBs seem to be a unique subclass of GRBs, and they are not the 
higher energy extension of the long GRBs. 



Key words: gamma rays: bursts — gamma rays: observations — methods: 
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1 INTRODUCTION 

Gamma-Ray Bursts (GRBs) are short and intense gamma-ray radiations from cosmo- 
logical distances. Much progress on GRBs and their afterglows has been made in the 
recent decade (Fishman & Meegan 1995; Piran 1999; van Paradijs et al. 2000; Cheng & 
Lu 2001; Meszaros 2002; Zhang & Meszaros 2004; Piran 2004). It is believed today that 
GRBs are produced by conical ejecta (jet) powered by central engines at cosmological 
distances. Numerous observations of this phenomenon have been successfully explained 
by the popular fireball models. However, this phenomenon is still much mysterious. A 
great number of puzzles are still to be solved. 

Unification of GRB population is one of the puzzles. Based on their distribution in 
the burst duration-harness ratio (HR) plane, Kouvcliotou et al. (1993) suggested two 
subclasses of GRBs, long GRBs and short GRBs separated at ~ 2 seconds (see also Qin 
et al. 2001). In recent years, a new subclass of GRBs, the X-ray flashes (XRFs), was 
also discovered (e.g., Heise et al. 2001). Whether or not the different subclasses of GRBs 
are the same phenomenon? How to provide a unified description for GRB population? 
These questions have been gotten much attention (Heise et al. 2001; Kippcn et al. 2003; 
Sakamoto et al. 2004a, 2004b; Lamb et al. 2003; Yamazaki ct al. 2003a, 2003b, 2004a, 
2004b; Lloyd-Ronning et al. 2004; Zhang et al. 2004; Dai & Zhang 2004; Liang & Dai 
2004). 

The peak energy of the vF v spectrum, E p , is an important quantity of GRB. Liang 
et al. (2002a) showed a limitation of E p with burst duration. This limitation cannot be 
explained by the current fireball model. It may represent a constraint on the fireball 
model. The narrowness of the E p distribution is also a great puzzle. The distribution 
based on a time-resolved spectral catalog of bright BATSE GRBs (Preece et al. 2000) 
shows that E p is narrowly clustered at 200-400 keV. It might be caused by the selection 
effect of the BATSE, which is sensitive to photons in energy band 50-300 keV. HETE-2 is 
sensitive to photons in an energy band of ~ 10—400 keV, and it is suitable for observations 
to XRFs. It is found that the spectra of XRFs are well fitted by the Band function (Band 
et al. 1993) with similar spectral indices as long GRBs (e.g., Kippcn ct al. 2003; Barraud 
et al. 2003), and they obey the same relations of E p to the equivalent-isotropic energy 
(Amati et al. 2002; Lloyd-Ronning & Ramirez-Ruiz 2002; Sakamoto et al. 2004a; Lamb et 
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al. 2003; Liang, Dai, & Wu 2004; Yonetoku et al. 2004) and to the jet energy (Ghirlanda, 
Ghiscllini & Lazzati 2004b; Dai, Liang, & Xu 2004). The duration distributions of long 
GRBs and XRFs are also similar (Heise et al. 2001). These facts suggest that long GRBs 
and XRFs are the same phenomenon, and XRFs are a lower peak-energy extension of 
long GRBs (e.g., Kippen et al. 2003; Sakamoto et al. 2004a). Therefore, XRFs extend 
the E p distribution to few keVs. Whether or not the E p distribution can be broadened 
to a higher energy is still uncertain (e.g., Piran 2004). Short GRBs tend to have harder 
spectra. Are short GRBs to be a higher energy extension of long GRBs? Ghirlanda, 
Ghiscllini & Celotti (2004a) presented an analysis to short GRBs. They found that the 
spectra of short bursts are well fitted by a single power law with an exponential cutoff 
at high energies. The statistics in the high energy channels of the spectra is too poor to 
constrain the high energy power law component of the Band function. They also found 
that spectral properties of short GRBs are similar to the first 1 second of long GRBs. 
Yamazaki et al. (2004b) proposed a unified description for long GRBs, short GRBs, and 
XRFs based on an off-axis jet model. If long and short GRBs are the same phenomenon, 
one may expect that short GRBs are a higher energy extension of long GRBs, and long 
GRBs, short GRBs, and XRFs form a sequence in the HR — E p plane, since the short 
GRBs tend to have a harder spectrum. In this work we investigate this issue with the 
current GRB samples observed by BATSE and HETE-2. 

This paper is arranged as follows. The method of analysis is presented in Section 2. 
The global spectral behaviors of the long GRB and XRF samples observed by HETE- 
2 are shown in Section 3, and those of the long and short GRB samples observed by 
BATSE, in Section 4. A discussion and conclusions are presented in Section 5. 



2 ANALYSIS METHOD 

It is well known that GRB spectrum can be well fitted by the Band function (Band et 
al. 1993), 



where A is a normalization parameter, Eq the break energy (in keV), a and (3 the low 
and high band spectral indices, respectively. For the case of (3 < — 2 and a > —2, the 
peak energy can be derived by E p = (2 + a)E , which corresponds to the energy at the 
maximum flux in vf u spectrum. For a spectrum with a > — 1 and (3 > — 2, the "true" E p 
lies at an unknown energy beyond the high end of the data, and only a lower boundary, 
i?b, to the energy of the high-energy power-law component characterized by (3 is known. 
The observed fluence at a given bandpass (E\, E2) as a function of E p can be given by 




E < (a — 0)Eo, 
E>(a-(3)E , 



(1) 



r E 2 



S Ei -e 2 (E p )= / EN(E)dE. 



(2) 
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The hardness ratio is defined as 



HR(E p ) ee 



(3) 



and its error caused by E p is (without consider the errors of a and (3 in our analysis), 



GRBs with similar spectral behaviors should trace out a sequence in the HR — E p 
plane. In this paper we investigate whether or not different subclasses of GRBs are in 
a HR — E p sequence characterized by a universal Band function by using the GRBs 
observed by BATSE and HETE-2. 

3 LONG GRBS VS. XRFS 

HETE-2 is sensitive to photons in the energy band of ~ 10 — 400 keV. It is suitable 
for observing GRBs and XRFs. We examine the sequence of long GRBs and XRFs in 
the HR - E p plane with the bursts observed by HETE-2. Of the 63 HETE-2 bursts 
(to the end of June, 2004), 34 bursts arc long GRBs and 22 bursts are XRFs. We in- 
clude only the long GRBs and the XRFs in our analysis. The values of E p , fluences 
(5 , 30-400koV and S2-30keVj m energy bands of 30 — 400 keV and 2 — 30 keV, respec- 
tively), and their errors for 40 bursts are taken from Sakamoto et al. (2004b). The 
other bursts are taken from HETE-2 burst home page 1 and G. Ricker's report (pri- 
vate communication), but no errors of E p and fluences are available for these. The errors 
of GRB030324 and GRB030723 given by Sakamoto et al. (2004b) are extremely large. 
We calculate the average percentages of the errors for the bursts from Sakamoto et al. 
(2004b) not includings GRB030324 and GRB030723, and obtain < a Ep /E p >~ 0.19, 
< CT5 3 o-4oo k cv/' S '30-40okoV >- 0.26, < cr S2 _ 30koV /5 2 - 3 okoV >- 0.09. Thus, for those bursts 
with no errors available, and for the two bursts, GRB030324 and GRB030723, we take 
<Je p = 0.19£p, cr S30 _ 400koV = 0.26S , 3O -4ookeV, $2-30keV - 0-095 2 -30k e v- The HR ob is cal- 
culated by HR ob = S , 3o-4ookov/5 , 2-30koV, and its error is given by Eq. (4) with measured 
errors of the fluences in the two energy bands. HR ob as a function of E p for the HETE-2 
bursts is shown in Fig. 1. From Fig. 1, we find that the GRBs and the XRFs seem to 
form a well-defined sequence. We fit this sequence with Eq. (3) by minimizing the x 2 > 




where 




E<(a- p)Eo, 
E > (a-/3)E . 



(5) 



* 2 = E 




(6) 



1 http://space.mit.edu/HETE/Bursts/ 
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Fig. 1 HR as a function of E p for the long GRBs (solid circles) and XRFs 
(open circles) observed by HETE-2. The solid curve is the best fit with Eq. (3) 
by fixing (3 = —5. 

where (T HR th is the error of HR caused by the uncertainty of E p , which is calculated by 
Eqs. (4) and (5). In our calculations, we find that this sequence is characterized by a; 
while f3 is numerically unstable (varying it from —5 to —10 gives comparable values of 
Xmin (~ H8)). So we fix /3 — —5, and then make a best fit to the sequence, and obtain 
a = —1.00 with Xmin/dof = 2-2. The best fit curve is also plotted in Fig. 1 (the solid 
curve). These results show that both the GRBs and the XRFs can be described by a 
universal Band function with E p ranging from few keVs to hundreds of keVs. The XRFs 
are at the lower end of the long GRB sequence, suggesting that XRFs are the lower 
energy extension of the long GRB sequence. 

4 LONG GRBS VS. SHORT GRBS 

In the present GRB sample, most of the short GRBs were observed by BATSE. To 
investigate whether or not both short and long GRBs can be characterized by a universal 
Band function, we use the long BATSE GRB sample (149 GRBs) given by Lloyd-Ronning 
& Ramirez-Ruiz (2002) and the short BATSE GRB sample (28 GRBs) presented by 
Ghirlanda et al. (2004a) . The fluences and their errors are available at the FLUX TABLE 
of the Current BATSE Catalog 2 . We calculate the hardness ratios by the observed 
fluences in channel 3 (110 — 320 keV) to those of channel 2 (55 — 110 keV), and their 
errors are derived by Eq. (4) with the errors of the fluences in channels 2 and 3. The HR 
versus E p plot for the long BATSE GRB sample is shown in Fig. 2. Since no E p errors of 
these bursts are available in Lloyd-Ronning & Ramirez- Ruiz (2002), we fit the long GRB 
sequence by Eq. (3) without considering the errors. The result shows that the minimum 
residual square, = Y,t( HR f ~ HRf 1 ) 2 , is 131 at (a, /?)=(-0.3, -2.05), which is also 
shown in Fig. 2. 
2 http:/ /cossc. gsfc.nasa.gov/batse/ 
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Log Ep 



Fig. 2 HR as a function of E p for the long BATSE GRBs. The solid curve is 
the best fit with Eq. (3). 




Log e p 

Fig. 3 HR as a function of E p for the short BATSE GRBs. The solid curve is 
the best with Eq. (3) with fixed (3 — —2.0 , and the dotted curve is the best fit 
for the long BATSE GRB sample. 



The HR versus E p plot for the short BATSE GRB sample is shown in Fig. 3. The 
spectrum of these short bursts is well fitted by a single power law with an exponential 
cutoff at high energies. Besides, the statistics in the high energy channels of the spectra is 
too poor to constrain the high energy power law component of the Band model (Ghirlanda 
et al. 2004a). We make a best fit to the short GRB sequence with Eq. (3). We find that 
f3 is ill-defined. So we fix (3 = —2 (similar to the long BATSE GRBs), and then make 
a best fit to the sequence, and derives a = —0.60 with X^m/dof — 1-1- The a value is 
quite similar to the mean of a in Ghirlanda et al. (2004a) (-0.58). The best fit curve is 
plotted in Fig. 3 (the solid curve). For comparison, the best fit for the long BATSE GRB 
sample is also plotted in Fig. 3 (dotted curve). It is found that the two best fits arc quite 
different. Most of the short GRBs are above the long GRB sequence, indicating that the 
short GRBs tend to have a larger HR than the long GRBs at a given E p . Short GRBs 
are unlikely to be a higher energy extension of long GRBs. 
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We have investigated the global spectral behaviors of long GRBs, short GRBs, and XRFs 
in the HR - E p plane with BATSE and HETE-2 observations. It is found that the long 
GRBs and the XRFs observed by HETE-2 seem to follow the same sequence in the 
HR — E p plane, with the XRFs at the low end of this sequence. We fit the sequence by 
a universal Band function, and find that this sequence is mainly fixed by a, while (3 is 
numerically unstable, with Xmin remaining at ~ 118 for (3 varying form —5 to —10. So 
we fix (5 — —5, and then make a best fit to the sequence, and derive a = —1.00 with 
Xmin/dof = 2.2. The sequences of the long and short GRBs observed by BATSE in the 
HR — E p plane are significantly different from that of the long GRBs. In the HR — E p 
plane, most of the short GRBs are above the long GRB sequence, indicating that the 
short GRBs tend to have a larger HR than the long GRBs at a given E p . A best-fit gives 
a = -0.30 and (3 = -2.05 for the long BATSE GRBs. For the short GRBs, a best fit is 
obtained with a = —0.60, with Xmin = I- 1 (ft is fixed at —2.0 since it is again numerically 
unstable). The a value of the short GRB sequence is significantly greater than that of 
the long GRBs. These results suggest strongly that the global spectral behaviors of the 
long GRBs and the XRFs are similar, and the XRFs are the lower energy extension of 
the long GRBs, while the global spectral behaviors of the short GRBs are different from 
the long GRBs: the short GRBs seem to be a unique subclass of GRBs, and they are not 
the higher energy extension of the long GRBs. 

The unified description of long GRBs and XRFs has been widely discussed. Zhang 
et al. (2004) showed that current GRB/XRF prompt emission/afterglow data can be 
described by a quasi-universal Gaussian-like (or similar structure) structured jet with 
a typical opening angle of ~ 6° and with a standard jet energy of ~ 10 51 ergs. Based 
on HETE-2 observations, Lamb et al. (2003) proposed that the uniform jet model nicely 
describes the prompt emission data of GRBs/XRFs, but this model fails to account for the 
afterglow jet break time data of GRBs. Liang & Dai (2004) found a bimodal distribution 
of the observed E p of GRBs/XRFs and suggested that two-component jet model (Berger 
et al. 2003) can explain this distribution, proposed that the two-component jet seems to 
be universal for GRBs/XRFs. The results of this paper further confirm the view that 
long GRBs and XRFs are the same phenomenon. 

We show that short GRBs are not the higher energy extension of GRBs, although 
they tend to have a harder spectrum than long GRBs. Please note that this result is 
based on the BATSE observations. Whether or not the E p distribution can be broadened 
to a higher energy is still uncertain. This possibility cannot be ruled out by BATSE 
observations. 

The observed association of long GRBs with star formation regions (e.g., Djorgovski 
et al. 2001), and the possible supernova components in afterglow light curves (e.g., Bloom 
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et al. 1999; Rechart et al. 2001; Stanek et al. 2003; Hjorth ct al. 2003) indicate that the 
central engines of the long GRBs might be the collapses of supermassive stars to black 
holes (e.g., Wooseley 1993; Bcrgcr, Kulkarni, & Frail 2003). Our results suggest that short 
GRBs are not the same phenomenon as the long GRBs and the XRFs. So far no afterglows 
of short GRBs have been observed. It is not clear whether this is an observational artifact 
or a real feature. Furthermore, we do not know of any direct evidence relating to the origin 
of short GRBs. The results in this paper might hint that short and long GRBs might 
come from different progenitors. Liang et all. (2002b, c) made the same argument based 
on the differences of the variability time scales and the total fluences of long and short 
GRBs. 

From Fig. 1, it is found that the sequences in the HR — E p plane are mainly charac- 
terized by the value of a, indicating that the GRB spectra in our analysis are dominated 
by photons with energies lower than E p . In our calculations, we find that Xmin is sensi- 
tive to a, but not to (3. The (3 value significantly affects the low end of the sequence at 
logi?p/keV < 1. However, there are only four XRFs in our GRB sample in this range. 
Thus, we cannot well constrain the low tail of the GRB sequence (and then the value of 
/?). In fact, the (3 value is ill-defined. This indicates that the statistics in the high energy 
channels of the spectra is too poor to constrain the high energy power law component of 
the Band model (Ghirlanda et al. 2004a). 

It should be noted that our best fit to the GRB sequence in the HR — E p plane is 
related to the definition of HR. Different HR definition can lead to different best fit 
results for the same GRB sample. This is the reason why the best fits shown in Fig. 
1 and 2 for HETE-2 long GRB sample and for BATSE long GRB sample are quite 
different. For the BATSE GRB sample, HR is calculated by the fluences in Channel 3 
and 2(encrgy bands 110-320 kcV and 55-110 keV) while the HR of the HETE-2 GRB 
sample is calculated by the fluences in the energy bands 30-400 keV and in 2-30 keV. 
Therefore, the best-fits to the HETE-2 GRB sample and to the BATSE GRB sample 
are different. However, this difference does not affect the results of our analysis because 
we separately investigated the long GRB-XRF sequence and the long GRB-short GRB 
sequence by the two instruments. 
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